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Abstract 

We study the formation and the evolution of velocity distribution tails for systems 
with weak long-range interactions. In the thermal bath approximation, the evolution of 
the distribution function of a test particle is governed by a Fokker-Planck equation where 
the diffusion coefficient depends on the velocity. We extend the theory of Potapenko et al. 
[Phys. Rev. E, 56, 7159 (1997)] developed for power-law diffusion coefficients to the case 
of an arbitrary form of diffusion coefficient and friction force. We study how the structure 
and the progression of the front depend on the behavior of the diffusion coefficient and 
friction force for large velocities. Particular emphasis is given to the case where the velocity 
dependence of the diffusion coefficient is Gaussian. This situation arises in Fokker-Planck 
equations associated with one dimensional systems with long-range interactions such as 
the Hamiltonian Mean Field (HMF) model and in the kinetic theory of two-dimensional 
point vortices in hydrodynamics. We show that the progression of the front is extremely 
slow (logarithmic) in that case so that the convergence towards the equilibrium state is 
peculiar. Our general formalism can have applications for other physical systems such as 
optical lattices. 

1 Introduction 

The study of Fokker-Planck equations is an important problem in statistical mechanics and 
kinetic theory pQ. In the simplest models, the diffusion coefficient is constant. However, Fokker- 
Planck equations with a diffusion coefficient depending on the velocity of the particles have also 
been introduced in physics. These equations usually describe the relaxation of a "test particle" 
evolving in a bath of "field particles" at statistical equilibrium when the particles interact via 
weak long-range forces. In that case, the diffusion coefficient is a function of the velocity of 
the test particle. For example, in his Brownian theory of stellar dynamics, Chandrasekhar [2] 
describes the evolution of the velocity distribution of a star in a cluster by a Fokker-Planck 
equation involving a diffusion and a friction. The coefficients of diffusion and friction are related 
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to each other by an Einstein relation and the diffusion coefficient decreases as v ~ 3 for large 
velocities. These results are similar to those obtained in plasma physics for the Coulombian 
interaction j^j. By using an analogy with stellar dynamics, Chavanis jU E] describes the 
relaxation of a test vortex in a thermal bath of field vortices by a Fokker-Planck equation (in 
position space) involving a diffusion and a drift along the vorticity gradient. The coefficients of 
drift and diffusion are related to each other by a form of Einstein relation involving a negative 
temperature and the diffusion coefficient is inversely proportional to the local shear created by 
the vortex cloud. For a Gaussian distribution of field vortices, the diffusion coefficient of the 
test vortex decreases with the distance as r 2 e~ Xr [HE]- Similarly, for the Hamiltonian Mean 
Field (HMF) model, Bouchet Sz Dauxois [01 E2] and Chavanis et al. ^T] find that the velocity 
distribution of a test particle satisfies a Fokker-Planck equation with a diffusion coefficient 
decreasing as e~^ v I 2 for large velocities. More generally, using the theory developed by Landau, 
Lenard and Balescu in plasma physics J2j and implementing a thermal bath approximation, 
one can obtain a general Fokker-Planck equation involving an anisotropic diffusion coefficient 
depending on the velocity of the test particle. This Fokker-Planck equation is valid for systems 
with weak long-range potentials of interaction. The preceding kinetic equations can be recovered 
as particular cases of this general Fokker-Planck equation jjj. 

For Fokker-Planck equations with a variable diffusion coefficient, the relaxation towards the 
Boltzmann distribution is slowed down, especially if the diffusion coefficient decreases rapidly 
with the velocity. One consequence is that velocity correlation functions can decrease alge- 
braically rapidly with time (instead of exponentially) as investigated by Bouchet & Dauxois 
[TU] in relation with the HMF model 1 . They therefore explain the observed algebraic tails 
of the velocity correlations functions in terms of classical kinetic theory without invocating a 
notion of "generalized thermodynamics". Here, we consider the relaxation of the system to- 
wards equilibrium from another point of view. We focus on the distribution function f(v,t) 
and study the structure and the evolution of the front formed in the high velocity tail. Our 
study is based on the approach of Potapenko et al. jT3] who studied this problem in the case 
where the diffusion coefficient decreases algebraically with the velocity. In the present paper, 
we generalize this approach to an arbitrary form of diffusion coefficient and study how the front 
position Vf(t) and the front structure evolve with time. In the case of a Gaussian or exponential 
decay of the diffusion coefficient with the velocity, we find that the progression of the front is 
extremely slow (logarithmic). This can leads to a sort of "kinetic blocking" or, at least, a 
"slowing down" of the relaxation. A similar confining effect was noted in the case where the 
diffusion coefficient depends on the density HI] , but this situation (related to the theory of 
violent relaxation) is more difficult to analyze since the kinetic equation is then nonlinear. 

The paper is organized as follows. In Sec. 01 we discuss different systems with long-range 
interactions (stellar systems, Coulombian plasmas, two-dimensional vortices, HMF model,...) 
which are described in the thermal bath approximation by Fokker-Planck equations with a 
variable diffusion coefficient. We also consider systems like optical lattices described by Fokker- 
Planck equations with constant diffusion coefficient but variable friction coefficient. In Sec. H3 
we generalize the theory of Potapenko et al. [T^] for an arbitrary form of diffusion coefficient 
and friction force. We provide general equations characterizing the structure and the evolution 
of the front for large times. In Sec. EJ we address the validity of our approach. In Sec. 
we consider particular applications of our general formalism for physically motivated Fokker- 
Planck equations. Analytical results are compared with direct numerical simulations of the 
Fokker-Planck equation. Particular emphasis is given to the case where the diffusion coefficient 

-"-In that case, the Fokker-Planck equation with Gaussian diffusion coefficient and linear friction can be 
transformed into a Fokker-Planck equation with constant diffusion coefficient and logarithmic potential which 
is known to exhibit power-law correlations (see, in particular, Appendix B of jl.'ij and |14|). 
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decreases with the velocity like an exponential or a Gaussian distribution. This is the situa- 
tion relevant for one-dimensional systems like the HMF model and for two-dimensional point 
vortices. Finally, in Sec. we investigate a class of Fokker-Planck equations for which our 
approach is exact for all times. The Appendices provide technical details and extensions of our 
main results. 



2 Examples of Fokker-Planck equations with a variable 
diffusion coefficient and friction force 

We consider a Hamiltonian system of N particles interacting via a weak long-range binary 
potential u(\r — r'\). These particles can be stars in stellar clusters, electrons or ions in a 
plasma, point vortices in two-dimensional hydrodynamics, particles located on a ring in the 
HMF model ... We assume that the cluster is homogeneous and in a steady state characterized 
by a distribution function /o( v )- In general, fo will be the statistical equilibrium state (thermal 
bath) but in certain cases it can be a slowly evolving distribution function. We introduce a "test 
particle" and denote by P(v,t) its velocity distribution function. Due to the interaction with 
the "field particles" , the velocity distribution will change and the test particle will acquire the 
distribution of the bath /o( v ) f° r t — > +00 (see [7j for more details). The general Fokker-Planck 
equation describing the relaxation ( "thermalization" ) of the test particle in the bath of field 
particles can be written as jjj: 

(i) 

where «(k) is the Fourier transform of u(r) and 

r k- ^ 



e(k,cu) = 1 + (2n) d u(k) / f—dv, (2) 

is the dielectric function. This equation can be obtained from the Lenard-Balescu equation by 
fixing the distribution function of the field particles /(v 1; t) to its static value /o( v i) [HI Ej- This 
(thermal) bath approximation transforms an integro-differential equation (Lenard-Balescu) into 
a differential equation (Fokker-Planck). The Lenard-Balescu equation was introduced in plasma 
physics for the Coulombian potential but it can apply to other systems with weak long-range 
interactions. Equations (|T])-P|) can also be obtained from the general expression of the Fokker- 
Planck equation by explicitly calculating the coefficients of friction and diffusion (first and 
second moments of the velocity increments) using the Klimontovich approach The fact 
that Eq. (JTJ) is linear does not imply that the distribution P(v,t) is close to equilibrium. The 
test particle approach is different from considering a small perturbation of the Lenard-Balescu 
equation around equilibrium. In the first case, we describe the evolution of a single test particle 
(or an ensemble of test particles that do not interact among themselves) in a thermal bath while 
in the second case one would describe the evolution of all the particles (the system "as a whole" ) 
close to equilibrium. 

If we consider that the field particles are at statistical equilibrium with the Boltzmann 
distribution /o( v ) ~ e -$ mv / 2 ; an d if we neglect collective effects taking |e(k, k-v)| = 1 (Landau 
approximation) we can rewrite the general Fokker-Planck equation (JTJ in the form [7j: 



BP 1 d 



dt t R dxL 



dP 

(r'U)[— + 2Px v 



(3) 
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where we have set x = (/3m/2) 1 / 2 v. The diffusion coefficient is proportional to the tensor 



with k = k/fc, and the quantity 



t 
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&^ 2 (2n) 



(4) 



d P m 

; ,3 



+00 



k d u(k) 2 dk 



(5) 



provides an estimate of the inverse relaxation time of the test particle toward the distribution 
of the bath. Here v m = (d/ ' (3m) 1 / 2 denotes the r.m.s. velocity and p the spatial density. A more 
general expression of the diffusion coefficient can be obtained by taking into account collective 
effects Note, however, that for v —>■ 00 (which is a limit that we shall be particularly 
interested with in the sequel), the expression (jlj is asymptotically exact since |e(k, k ■ v)| — > 1 
for |v| — > +00. Note also that for weak long-range potentials of interaction for which our 
approach is valid, and in the Landau approximation, the precise form of the potential u(k) 
only determines the timescale of the relaxation, through Eq. not the form of the kinetic 
operator. Therefore, in the Landau approximation, the expression of the diffusion tensor as a 
function of the velocity only depends on the dimension of space d. 
In dimension d = 3, the diffusion tensor can be written 



G*™ = (Gii 



1 



G 



1 



G 1 5^, 



(6) 
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where Gu and G± are the diffusion coefficients in the directions parallel and perpendicular to 
the velocity of the test particle. They are explicitly given by 



with 



where 



G(x) 



2tt 3 / 2 



G(x), 



G 



2tt 3 / 2 



x 



erf (a;) — G(x)], 



(7) 



2 1 



7T X' 



t 2 e~ t2 dt 



2x 2 



erf (x) 



erf(x) 



~ t2 dt, 



2x _ T 2 

"K 



(9) 



is the error function. If we consider spherically symmetric distributions, noting that dPj dx^ 
(l/x)(dP/dx)x fM and G^x" = Gnx^ we obtain 



dP 
~dt 



1 1 d 

t R x 2 dx 



( dP 

x 2 G\\(x)[— + 2Px 



(10) 



For the gravitational potential, this Fokker-Planck equation has been studied by Chandrasekhar 
in his Brownian theory of stellar dynamics 0. It has also been considered in plasma physics 
as an approximation of the Landau equation valid for sufficiently large times [3 E3- We note 
in particular that the diffusion coefficient G\\{x) decreases algebraically like x~ 3 for x — > +00. 

Alternatively, if we consider one dimensional systems (d — 1), the general Fokker-Planck 
equation (JTJ) simplifies into jjj: 



dP 
~dt 



d_ 

dv 
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where D(v) is given by 



D(v) = An 2 mf (v) f °° dk . (12) 
Jo \e{k,kv)\ 2 

We note that the distribution function P(v,t) of the test particle relaxes towards the distribu- 
tion of the bath fo(v) on a timescale of the order Nt^, where to is the dynamical time [JJ. If 
we neglect collective effects, or consider the limit of large velocities, we find that the diffusion 
coefficient is given by 

r+oo 

D(v) = 4ir 2 mf (v) / dkkuik) 2 . (13) 
Jo 

It is proportional to the distribution function of the bath fo(v). In particular, if the field 
particles are at statistical equilibrium with a Gaussian distribution, the diffusion coefficient 
decreases like e~ l3mv I 2 . This type of Fokker-Planck equations apply for example to the HMF 
model P|llO|ITT] which can be viewed as the one Fourier component of a one-dimensional plasma 
(or self-gravitating system) [TT]. More generally, these Fokker-Planck equations ([lip are valid 
for a wide class of one dimensional systems with long range interactions [JJ. We note that for 
one-dimensional systems the Lenard-Balescu collision term cancels out so that the distribution 
function of the field particles f(v,t) does not evolve, i.e. df/dt = 0, on a timescale of order 
Ntp- Since, on the other hand, the relaxation time of a test particle toward the distribution 
of the bath is of order Nto, this implies that we can assume that the distribution of the field 
particles is stationary f(v,t) = fo(v) when we study the relaxation of a test particle. This 
is true for any distribution function fo(v) that is a stable stationary solution of the Vlasov 
equation [TUJ HU [7j- This is not true in higher dimensions d = 2 and d = 3, except for the 
Maxwellian distribution / e (v), since the distribution of the field particles /(v, t) changes on a 
time Nto as it relaxes towards the statistical equilibrium state / e (v). Furthermore, even if we 
assume that the distribution of the bath /o(v) is approximately stationary, the distribution of 
the test particle P(v, t) will not relax towards /o( v ) f° r large times except if /o( v ) is Maxwellian 

m 

A Fokker-Planck equation with a space dependent diffusion coefficient has been introduced 
by Chavanis in [3J E] to describe the relaxation of a test vortex in a "sea" of field vortices 
with vorticity profile u>o(r). For an axisymmetric distribution, the Fokker-Planck equation for 
P(r,t) can be written 0[7j : 



dP _ 1 d 

dt r dr 



(14) 



with a diffusion coefficient 



D{r)= 8^H lnAMr) ' (15) 

where S(r) = rQ' (r) is the local shear created by the field vortices (£lo(r) represents the angular 
velocity related to the vorticity by ui (r) = (l/r)(r 2 fi )')- F° r a vorticity profile u (r) = Ae~ Xr 
of the field vortices, it is easy to see that the diffusion coefficient of the test vortex decreases 
like D(r) ~ r 2 e~ Xr2 for r -> +oo EJ. 

Finally, Fokker-Planck equations with diffusion and friction coefficients depending on the 
velocity can occur in many areas of physics. For example, the motion of atoms in a one- 
dimensional optical lattice formed by two counter propagating laser beams with linear perpen- 
dicular polarization can be described, after spatial averaging, by a Fokker-Planck of the form 
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dW _ d 
dt dp 



dW 

D(p)—-WK(p) 



with 



K{p) 



ap 



1 + (p/p c 



\2> 



D{p) = D c 



1 + (p/p c 



(16) 



(17) 



For p — > +oo, D(p) — > Dq and K(p) ~ —ap 2 c /p. This corresponds to a logarithmic potential 
U{p) ~ (otp 2 c / D )\np defined by K{p)/D{p) = —U'(p). Equation (fTo^) belongs to the general 
class of Fokker-Planck equations that we shall study in the sequel. The case of a logarithmic 
potential is treated specifically in Sec. 15 .HI 



3 General solution of the problem 

The various examples discussed previously prompt us to study Fokker-Planck equations with a 
diffusion coefficient and friction force depending on the velocity. In particular, we can wonder 
how the distribution function f(v,t) approaches the equilibrium distribution. This problem 
has been investigated by Potapenko et al. ^S] in the case of 3D plasmas where the diffusion 
coefficient is a power law. These authors found that the asymptotic behavior of the velocity 
distribution tail has a propagating wave appearance. The high velocity tail develops a front 
at which the distribution function drops to zero. This front Vf(t) progresses with time and 
goes to Vf(t) — > +oo for t — > +oo. The profile of the front also deforms itself as time goes 
on. However, if we use an appropriate system of coordinates, it can be expressed in terms of 
the error function. We shall here formulate the problem in a general setting, for an arbitrary 
form of diffusion coefficient and friction force, and we shall investigate how the evolution of the 
front and the profile of the high velocity tail distribution depend on the form of the diffusion 
coefficient. 

Let us consider the Fokker-Planck equation 



df 1 d 



dt v d 1 dv 

where D(v) > and U(v) are arbitrary functions. If the following zero flux condition 



v d ~ l D{v) \M + fU'{v) \ - (19) 

when v — > oo is fulfilled, then the stationary solutions of this Fokker-Planck equation take the 
form 

f e (v) = Ae- u ^, (20) 

where A is a constant (normalization). The Fokker-Planck equation (|18J) decreases the free 
energy F = E — S where E = J fU(v)dv is the energy and S = — J fin fdv is the Boltzmann 
entropy (the temperature has been included in the potential U). Indeed, one has 

D f s f + f s JL) 2 dv < , (21) 



/ \dv dv J 
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Therefore, if F is bounded from below, the distribution will converge towards the equilibrium 
state ()20|) for t — > +00. We want to analyze the propagation of the front in the high velocity 
tail of the distribution function. Thus, we set 



f(v,t) = f e (v)u(v,t). 



(22) 



For sufficiently large times, the core of the distribution function will have reached its asymptotic 
value ()20j) so that u ~ 1 in that region. On the other hand, for sufficiently large velocities, 
the distribution has not relaxed yet and u = 0. Therefore, we expect the formation of a front 
at a typical velocity value ~ Vf{t) where the function u(v,t) passes from u — 1 to u — 0. On 
this phenomenological basis, u(v,t) is the relevant function to consider in our "travelling front" 
analysis. Its exact evolution is governed by 



du 
~dt 



l_d_ 



,d-l 



D(v)^)-D(v)U'(v) 



du 
dv ' 



(23) 



which is obtained from Eqs. (jl8)) - ([22|) . If we perform the change of variables dx/ dv = 1/ \J D{v). 
we obtain the equivalent equation 



d 2 u 
dx 2 



du 
~dt 



1 D'(v) d 
' + — 



-^D{v)-U'{v)^/D{v) 



du 
dx ' 



(24) 



where v = v(x) must be viewed as an implicit function of x. If we introduce the velocity field 



V(v) = 

or, equivalent ly, 

V{v) = 

Eq. (j24*|) can be re-written as 



'DM 



U'(v) 



d-1 1 



-(InD)'(v) 



(25) 



d 



du .du 

m +v{v) d- x 



-l e -U(v) D l/2, 



d 2 u 



(26) 



dx 2 



(27) 



The structure of this equation is clear. The right hand side corresponds to a diffusion and 
the left hand side corresponds to an advection by a velocity field V(v). This velocity field (in 
velocity space) governs the evolution of the front. To get a physical insight into the problem, 
let us first neglect the diffusion term. The resulting equation 



du ,du 

Tt + V(V) T* 







can be solved with the method of characteristics. Writing 

1 dv 



dx y. . dx dv 
dt dv dt 



r D{v)dt 



(28) 



(29) 



we get 



^ = ^W)V(v) = D{v) 



U'(v) 



-(InD)'(v) 



(30) 
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This equation determines the evolution of the front Vf(t). In the extreme approximation where 
the diffusion is neglected, the profile of the front is given by a step function u(x, t) = r)(x—Xf(t)) 
where rj is the Heaviside function (i](x) = 1 for x < and r](x) = for x > 0). As we shall see, 
the diffusion term will smooth out this profile. 

To take into account the effect of diffusion, we return to Eq. (|27|) and perform the change 
of variables 

z = x — Xf(t), u(x,t) = <f>(z,t), (31) 
where the function Xf(t) is defined by Eq. (|29|). Substituting this in Eq. (J27J), we obtain 

dd> d 2 (k ;! • 



[V(v)-V(v f )}£. (32) 



dt dz 2 dz 



So far, no approximation has been made so that Eq. (}32jl is exact and bears the same infor- 
mation as the initial Fokker-Planck equation (fTHj) . It is just written in a more convenient form 
which will allow us to examine the situation in the region of the front. Indeed, far from the 
front the profile is stationary (for sufficiently long time) and Eq. (1321) is automatically satisfied 
as = 1 for v <C Vf(t) and = for v ^> Vf{t). If we consider values of the velocity that 
are close to f/(t), we can expand the term in brackets in Taylor series (the validity of this 
approximation will be studied in Sec. ^). Keeping only the first term in this expansion, we get 

d(j) d 2 (j) d(j) 

m = d^- 9it)z ^ (33) 

where 

g(t) = V'(v f (t))^D(v f (t)). (34) 

For future convenience, we set r = 2t and define h{r) = g{r/2). Therefore, the forgoing 
equation becomes 

l '"*- HT) M\ (35) 



dr 2\ dz 2 dz t 

The general solution of this equation, for an arbitrary initial condition, is given in Appendix 
lAl Here, we look for particular solutions of the form <P(z,t) = $(z/x(r)). In Appendix IA"| we 
derive the condition under which such solutions describe the asymptotic long time behavior of 
the system (independently of the initial condition) and we check that this condition is fulfilled 
for all the explicit examples that we shall investigate in the following. To describe more general 
situations, one must use the results of Appendix IA1 Inserting the ansatz <p(z,r) = $(z/x(r)) 
in Eq. (j33|) leads to 

$" + (2 X X - h X 2 )x& = 0. (36) 

The variables of velocity and time separate provided that the term in parenthesis is a constant 
that we can arbitrarily set equal to 

2 X X - h X 2 = 2. (37) 

Then, $ is the function 

$(*) = — / e~y 2 dy, (38) 

7T L. 



8 



connected to the error function by = |(1 — erf (a;)). It satisfies $(— oo) = 1, $(0) = 1/2 
and <3>(+oo) = so it reproduces the expected properties of the front (it can be seen as a 
smooth step function). On the other hand, solving for \ 2 m Eq. ()37|) and taking x(l) = 0, we 
finally obtain 



e [H(T)-H(T>)] dT , 



(39) 



where H is a primitive of h with H(l) = 0. The function (p(z, r) = $(z/x(t)) with (jHHj) and 
f!39|) is the solution of Eq. (J33j) for all times, corresponding to the initial condition <f)(z, 1) = r](z) 
where rj is a step function. It also governs the long time behavior of the system for a large class 
of initial conditions (Appendix |XJ) . 

In particular, for h(r) = 7/r, we have 



1 ( d 2 (f) 7 dcj) 
dz 2 r dz 



dr 2 

and we recover the solution given in [E], namely 



(40) 



Z,T 



' z ( 1 -7 



1/2- 



For 7 



V2 V r - r7 , 

•1, a value that will frequently occur in the following examples, we have 

z 



[z,t) = $ 



-1/2 



for r > 1. 



(41) 



(42) 



In terms of the function g, these analytical solutions correspond to g(t) = X/t. We come back 
to the original variables by setting r = It and 7 = 2A. 



4 Condition of validity 

We shall now investigate in greater detail the ability of our approach to describe the front 
structure. First, we note that, within our approximations, 



u(v f (t),t) = ${0) = -, 
so that Vf(t) gives the position of the half- height profile. On the other hand, noting that 



(43) 



V(v) - V(v f ) = V'(v f )(v - v f ) + -V"(v f )(v - v f f + 



(44) 



for v — > v f, our approximation ()33|) will be valid in the range of velocities where we can neglect 
the second term in the Taylor expansion with respect of the first. This corresponds to 



\v — Vf \ <C 



2V(v f ) 



V"(v f ) 



(45) 



In this range of velocities, our approach is accurate. Now, it will provide a good description of 
the whole front if this range \v — Vf \ is larger than the typical front half-width Af(t)/2. This 
condition can be written 



A f {t) « 



W (v f ) 



(46) 
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Now, the front width can be estimated by 

A/(t) 



du 
dv 



(v f (t),t) 



-i 



(47) 



Within our approximation, this can be finally rewritten as 



A / (t) = v ^x(2t) V / ^Mt)). (4 
Therefore, our approach will provide a good description of the front structure if 



e(t) = 



4 V'(v f )/V»(v f ) 



< 1. 



(49) 



Since we shall be interested by the large time limit, it is particularly important to know the 
asymptotic behavior of the function e(t) for t — > +oo. This has to be considered case by case 
(see Sec. EJ). 

The other assumption made in our study is that the long time behavior of the system is 
described by 4>(z,t) = $(z/x(r)) with and (jHSJ)- In Appendix \K\ we show that this is the 
case if d>(z, 1) — > 1 for z — ■> — oo, <b(z, 1) — > for z — > +oo and if 



H 9 (t) 



for r — > +oo. In particular, for h{r) = 7/r, we have 

H(t) =7 In r, 

and 



X 2 (r) 



r — t 1 



H2{r) 



-1-7 



1 — 7 1 — 7 

so that the criterion (|5()jl is met if 7 < 1 and not met if 7 > 1. 



(50) 



(51) 



(52) 



5 Particular examples 

We shall now discuss explicitly some particular examples of physical interest and compare our 
theoretical predictions with direct numerical simulations of the Fokker-Planck equation. 

5.1 Quadratic potential 

We first consider the case of a quadratic potential U (v) = v 2 /2 leading to the Kramers equation 

df 1 d 



dt v d ~ l dv 



(53) 



with a variable diffusion coefficient D(v). If we assume a zero flux condition (|19|h then the 
stationary solutions are the Maxwellian distributions 

f e (v) = Ae-$. (54) 

We now consider various forms of diffusion coefficient and use the general theory developed 
in Sec. H3 to characterize the front profile and its evolution. We shall only give asymptotic 
expressions which are valid for large velocities and large time. This will be implicit in all the 
following calculations. 
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5.1.1 Power-law decay 

For the diffusion coefficient 

we get 



D(v) ~ v~ 



x ~ 



and 



a + 2 

git) 

9(t) 



a + 2 

-V 2 . 



a > 0, 



Vf(t) ~ (at) 



l/a 



a 1 



2a f 
d-21 



if a^2 
if a = 2. 



(55) 



(56) 



(57) 



2 t 2! 

For a = 2, we have assumed that the subdominant corrections to the D(v) ~ t>~ 2 behavior 
are of order v ~ 4 or smaller. Note that if a = 2 and d = 2, then g = and the Fokker-Planck 
equation can be solved exactly (see Sec. EJ). 

For a ^ 2, we find that h = j/t with j = 2/a — 1. The front profile is given by 



u(v,t) ~ $ 



a + 2 
V 2 



(at) 



a + 2 
2a 



1-7 



w-i 



1/2- 



(5f 



(a + 2)^/2 V 1 "^)" 

The criterion (|5()jl is fulfilled only for a > 1, so that the above function provides the correct 
asymptotic behavior of the solution, for any initial condition, only in that case. For a < 1, it 
can however provide the correct asymptotic behavior if the initial condition <f>(z, 1) is a step 
function (see Appendix |XJ). Equation (j3Hj) returns the results obtained in [T^]. Of course this 
formula is written for a / 1, but it also provides the expression of the solution for a = 1 by 
passing to the limit a — > 1 yielding: 



u(v,t) ~ $ 



a 3 

t>2 — f2 

3(tln(2t))V2 



(59) 



On the other hand, for a = 2 we find that /i = 2(d — 2)/r 2 yielding x 2 (r) ~ 2r and H 2 (t) ~ r 
for r — > +oo. We find that the criterion (|5U|1 is satisfied and that the front profile is given by 

• 2 2t 



u(v,t) ~ $ 



V 



4t l/2 



(60) 



which turns out to be consistent with Eq. (|3S)l. 

Concerning the validity of this approach in relation with the criterion (|49|1. we have 



e(t) 

e(t) ~ — 



e(t) 
e(t) 



if a > 1 and a^2, 

if a < 1, 

T(lnt)V 2 

- if a = 1, 

if a = 2 (d ^ 2), 



(61) 



4 

3v^F 



-1/2 



where A a = ( v / 2¥/8)a 1 - 1 / Q (a - l)' 1 ' 2 and fi Q = ( v ^F/8)a 1 - 1 / a (l - a)- 1 / 2 2 1 / Q - 1 / 2 . Therefore, 
the condition e <C 1 is always fulfilled for sufficiently large times. For sake of illustration, we 
show the case a = 3 for d = 3 (plasmas and stellar systems) in Figs. [Hand El These results are 
obtained by solving numerically the Fokker-Planck equation (|53|) starting from a step function: 
f(v, t = 0) = 3/47T if v < 1 and f(v , t = 0)=0iff>l (water-bag). In the simulation we have 
adopted the expression (J2J)-a of the diffusion coefficient which reduces to Eq. (fHHJ) with a = 3 
for v — ► +oo. Other examples are given in [To] . 
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Figure 1: Evolution of the front profile u(v, t) for a power-law diffusion coefficient with a = 3. 
For sufficiently large times, we get a perfect agreement with the theoretical profile Here 
and in the following, the straight lines correspond to the numerical simulation and the dashed 
lines to the theoretical prediction. 




Figure 2: Evolution of the front position f/(t) defined by Eq. (J43)) for a power-law diffusion 
coefficient with a = 3. For sufficiently large times, it coincides with the theoretical prediction 
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5.1.2 Gaussian decay 

We now consider a diffusion coefficient of the form 

D(v) ~ e"^ 2 . (62) 
Considering the evolution of the high velocity tail, we have 

x~ e% y2 dy. (63) 



o 



For large v, we obtain the relation (see Appendix |B|) 



xrv—efr*. (64) 
The position of the front Vf(t) is determined by 

— dy ~ 2( 7 + l)t. (65) 

The integral can be expressed in terms of the exponential integral Ei(x). For large times, we 
get 

—= -~2( T +l)t. (66) 

To leading order, we have 

An A 1/2 , , 

Vf(t) ~ (^J • ( 6? ) 

We note that the evolution of the front is extremely slow (logarithmic). On the other hand, we 
find that 

g(t) ~ ~, (68) 

implying that the criterion (J5U)) is fulfilled. In order to determine the front profile, we need 
first to evaluate Xf(t). An equivalent for t — > +oo is obtained by combining Eqs. (jMj) and 
leading to 



*/(«) - t/^- (») 



Therefore, the front profile is given by 



u(v,t) ~ $ 



ri e ]» 2 - ,/ 2(7+1)* n 

7« y 7 



(70) 



'2t 

Concerning the validity of this approach in relation with the criterion flU} , we find that e(t) — > 
i(t?i) 1//2 ^ or ^ ~~ > ~^ oc so ^^ a ^ our approximations are marginally valid: e(t) does not go to zero 
but it does not diverge with time neither. 



13 




Figure 3: Short time evolution of the distribution function f(v,t) for a Gaussian diffusion 
coefficient. This figure shows the slow depletion of the high velocity tail due to the rapid 
decrease of the diffusion coefficient for v Vth — 1. 



We have performed numerical simulations of the Fokker-Planck equation ()53|) with the 
diffusion coefficient (J52*|l with 7 = 1/2. As discussed in Sec. EJ this equation describes the 
evolution of a "test particle" in a bath of "field particles" at statistical equilibrium (with 
Maxwellian distribution function /q ~ e~~ v I 2 ) for one- dimensional systems (d = 1) with long- 
range interactions such as the HMF model. In Fig. El we show the "short" time evolution of 
the distribution function f(v,t) starting from a step function: f(v, 0) = 1/3 for v < 3 whose 
width is larger than the thermal speed v t h = 1 (dashed line). Since D{y ) rapidly decreases with 
the velocity for « > % = 1, the high velocity component of the initial condition takes time 
to be depleted. Indeed, the core of the distribution rapidly reaches a Maxwellian distribution 
while the tail keeps the memory of the initial condition and remains flat for relatively long 
times. For the parameters of the simulation, this can even create a non-monotonic distribution 
for short times as shown in Fig. E] (see Appendix IU|) . 

In Fig. H]we show the evolution of the normalized distribution function u(v , t) for short and 
large times. In that case, we start from an initial condition: f(v , 0) = 1 if v < 1 and f(v , 0) = 
if v > 1. In Fig. El we focus on the evolution of the front for large times and we compare 
the result of the numerical simulation with the theoretical prediction. In Fig. we compare 
the evolution of the front displacement Vf(t) with the theoretical prediction. The numerical 
results are in fair agreement with the theory although its domain of validity was expected to 
be marginal in that case according to our estimates. This fair agreement may be explained by 
the relatively small value of e(+oo) = ^(7r/3) 1//2 ~ 0.255... for 7 = 1/2 even if this parameter 
does not strictly tends to zero for t — > +00. We note that for t — > +00 and v <C Vf(t), we have 
u ~ $(— sj (7 + l)/7) = $(— v / 7r/4e) which has not converged to u = 1 precisely because e ^ 0. 
However, for e(+oo) ~ 0.255..., we find that u ~ 0.993 which is very close to one. 

In relation with the HMF model, we would like to recall that the present approach describes 
the relaxation of a given test particle immersed in a bath of field particles at statistical equi- 
librium. The relaxation is due to finite iV effects (correlations). This does not describe the 
evolution of the distribution function of the system as a whole. The clearest reason is that 
the Fokker-Planck equation resulting from a thermal bath approximation does not conserve 
energy contrary to the microcanonical evolution of a Hamiltonian system. In the collisionless 
regime, the distribution function f(6,v,t) of the HMF model is governed by the Vlasov equa- 
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Figure 4: Evolution of the normalized distribution function u(v,t) for a Gaussian diffusion 
coefficient for short and large times. For short times the normalized distribution function forms 
a bump which slowly disappears as the core of the distribution becomes Maxwellian. For large 
times the function u(v,t) has a front structure. 




Figure 5: Evolution of the front profile u(v,t) for a Gaussian diffusion coefficient. For suf- 
ficiently large times, we get a fair agreement with the theoretical profile ([70)1 although the 
validity of our approach was expected to be marginal in that case. 
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Figure 6: Evolution of the front position Vf(t) defined by Eq. (jjHj) for a Gaussian diffusion 
coefficient. For sufficiently large times, we get a fair agreement with the theoretical prediction 

(jnnj). 

tion coupled to the mean- field potential $(0,t) = -± J cos(t9 - 9')f(9',v',t)d9'dv' produced 
by the particles. The Vlasov equation can experience a process of violent relaxation and con- 
verge toward a Quasi Stationary State (QSS) on the coarse-grained scale [IB]. Some dynamical 
theories of violent relaxation JE] based on a Maximum Entropy Production Principle (MEPP) 
propose to model the evolution of the coarse-grained distribution function by a generalized 
Fokker-Planck equation of the form (for a water-bag initial condition): 



df df „fif d f ,„ 



^- + P(t)f( Vo -f)v 



(71) 



where f3(t) evolves so as to conserve energy (see ^HJ for more details). An important point is 
that the diffusion coefficient is not constant but is related to the correlations of the fine-grained 
fluctuations. It can depend on position, velocity and time and can be very small in certain 
regions of phase-space and for large times. The vanishing or smallness of the diffusion coefficient 
can slow down the dynamics and lead to a confinement of the distribution function in phase 
space which may be qualitatively similar to what is shown in Fig. El However, the dynamical 
equation (|7T|) and its physical interpretation are different from Eq. (fTTj) . so that their similarity 
is, at most, an analogy. 



5.1.3 Exponential decay 

For the diffusion coefficient 

we get 

The evolution of the front is given by 



D(v) ~ e 



2 2, 

x ~ — e 2 

7 



e jv f 



t. 



(72) 



(73) 



(74) 
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Figure 7: Evolution of the front profile u(v,t) for an exponential diffusion coefficient. For 
sufficiently large times, we get an excellent agreement with the theoretical profile (fT%|) . 



To leading order, we have 



Vf(t) ~ — hit, 



(75) 



(76) 



which shows that the progression of the front is again very slow. We also have 

m ~ -I, 

implying that the criterion (|5()jl is fulfilled. In order to determine the front profile, we need 
first to evaluate Xf(t). An equivalent for t — > +oo is obtained by combining Eqs. (|7Sj) and ()74|) 
leading to 



^--(tlnt)^. 



7 



(77) 



However, since the evolution with time is slow, we shall work with the more precise expression 
Xf(t) ~ ^(tln(t hit)) 1 / 2 obtained by keeping the term of next order. With this expression, we 
find that the front profile is given by 



u(v,t) ~ $ 



V2 



e* v - (t\n(t\nt)y/ 2 



-ft 1 / 2 



(7f 



Concerning the validity of this approach in relation with the criterion ()49j). we find that e(t) ~ 
|v / 27r(lnt) -1 ^ 2 for t — > +oo so that our approximations are valid for sufficiently large times. 
Note that the decay is slow with time (logarithmic) but the prefactor |v^27r = 0.313... is 
relatively small (for 7 = 1) which can explain the very good agreement with the numerics even 
for moderate timescales. Figures [7| and |H] are obtained by solving the Fokker-Planck equation 
(153)) with the diffusion coefficient (|72j) with 7=1, starting from an initial condition: f(v , 0) = 1 
if v < 1 and / = if v > 1. 



5.1.4 Stretched exponential 

For the diffusion coefficient 



D(v) ~ e 



5>0, 



(79) 
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Figure 8: Evolution of the front position Vf(t) defined by Eq. P3|) for an exponential diffusion 
coefficient. For sufficiently large times, we get an excellent agreement with the theoretical 
prediction (|71jl . 



we get 



x 



ez yS dy. 



(80) 



We shall briefly discuss how the results depend on the value of 5. A more detailed analysis 
can be carried out as in the previous sections where the cases 5 = 1 and 5 = 2 are explicitly 
considered. 

For 5 < 2, the front evolution is given by 



To leading order, we get 



St. 



ln/V 

7 



On the other hand, 



9{t) ~ ~ 



(82) 



J3) 



implying that the criterion (jHOj) is fulfilled. Concerning the validity of this approach in relation 
with the criterion (|49jl . we find that 



e(t) 



V2n5 flat 



7 



-(2-S)/25 



(84) 



so that the criterion is satisfied for sufficiently large times (but slowly). 
For 5 > 2, the front evolution is given by 



e< v f 



( 7U 5)2(«5-l)/5 2 7 



(85) 
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To leading order, we get 

1/5 



(86) 



On the other hand, 



g(t) ~ ~, (87) 

implying that the criterion (}50|) is fulfilled. Concerning the validity of this approach in relation 
with the criterion (|49p. we find that e — > ^ so that this approach is marginally valid. 

5.2 Linear potential 

We now consider the case of a linear potential U (v ) = 'jv leading to a Fokker-Planck equation 
of the form 



df 1 d 



dt v d ~ l dv 
The stationary solution is 

f e (v)=Ae~^. (89) 
For a diffusion coefficient decreasing algebraically with the velocity 

D(v) ~ v~ a , a > 0, (90) 

we get 

2 a+2 



X 



a + 2 
and 



V s *-, v f (t) ~ [(a + l) 7 t] 1/(a+1) , (91) 



^(t) I (92) 

yw 2(a + l)t' V ; 

implying that the criterion (|5Uj) is fulfilled. We also find that e(i) oc t -1 /^^ 1 )) so that the 
validity criterion of our approach is always fulfilled for sufficiently large times. 

5.3 Logarithmic potential 

In this subsection, we consider the case of a constant diffusion coefficient D(v) = 1 and a 
logarithmic potential U(v) = (a/2) ln(l + v 2 ) in d = 1 leading to a Fokker-Planck equation of 
the form 

This type of Fokker-Planck equations arises in the study of optical lattices Ej. The 
stationary solution of this equation is of the form 

AM = (TT^' (94) 
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Figure 9: Evolution of the front position Vf(t) for a logarithmic potential with a = 10. We get 
a fair agremment with the Lambert function. 



This is similar to a Tsallis distribution with q = (a — 2) /a [19 . However, it arises here from 
a linear Fokker-Planck equation (associated with the Boltzmann entropy) with a logarithmic 
potential, instead of a nonlinear Fokker-Planck equation (associated with the Tsallis entropy) 
with a quadratic potential [2H1 ED 122] • The distribution ((OH) is normalizable provided that 
a > a crit = 1. The case a = 2 corresponds to the Lorentzian. The evolution of the front Vf(i) 
satisfies 

dt ~ i+ v y [Jb) 

The exact solution is given by vje v f = e 2at+c where C is a constant of integration. It can 
be written v 2 (t) = W(e 2at+C ) where W(x) is the Lambert function which is solution of the 
transcendental equation W / (x)exp[W / (x)] = x. For t — ► +oo, we get Vf(t) ~ {2at) l l 2 and 
g(t) ~ —l/(2t). Therefore, the front profile is given by 



for large times. Concerning the validity of our approach with respect to the criterion (|49jl. we 
find that e -> (^) x / 2 for t -> + oo so that our approach is marginally valid. The theoretical 
prediction is not very good for a = 2 but the agreement improves for large values of a ^> 
a cr it = 1 for which e is reduced. In particular, the results of numerical simulations performed 
with a = 10 are shown in Figs. El and Concerning the evolution of the front, we find a 
good agreement with the Lambert function except that the measured exponent in Fig. El is 
~ 20.6 instead of 2a = 20. Therefore, the front increases a little bit faster than the theoretical 
prediction. This is confirmed in Fig. El which shows the evolution of the front profile. We 
note, however, the relatively good agreement between direct numerical simulation and theory. 
The slight discrepency is due to the finite value of e ~ 0.28. 
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Figure 10: Evolution of the front profile u(v,t) for a logarithmic potential with a = 10. 

5.4 Fokker-Planck equations associated with one-dimensional sys- 
tems with long-range interactions 

We now consider the case of Fokker-Planck equations describing one dimensional systems with 
long-range interactions. In that case, the diffusion coefficient D ~ fo{v) for v — > +oo and the 
potential U{v) = — \nf (v) are expressed in terms of the distribution of the bath fo(v). The 
corresponding Fokker-Planck equation can be rewritten as (see Sec. EJ): 



dt dv V dv dv 



The general relations obtained in Sec. El take the simpler forms 
dx 1 3 ft(v) 

Tv ~ v{v) ~ 



2 VM 



dv f 3 
■W ~ -2 foiVf) - 



(97) 



(98) 



The case where fo(v) is a Gaussian distribution (thermal bath) has already been studied in 
Sec. 15.1.21 We shall consider other examples where fo(v) is not necessarily the statistical 
equilibrium state. This description still makes sense if /q is a stable stationary solution of the 
Vlasov equation because the relaxation of the system as a whole is slower than the relaxation 
of a test particle in a fixed distribution. 

5.4.1 Exponential distribution 

For the exponential distribution 



fo ~ e 



we obtain 



and 



x ~ — e 2 

7 



i , a 



Vf ~ — In -7 i 



(99) 



(100) 



g(t) 



2t' 



lor 
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implying that the criterion (J5U|) is fulfilled. The front profile is given by 



u(v,t) ~ $ 



y/2 



e/2 



v _ f' 3 2^1/2 



-ft 1 / 2 



(102) 



Concerning the validity of this approach in relation with the criterion (|49p. we find that e — > 
i(l) 1 / 2 ~ 0.256... for t — > +oo so that our approximations are marginally valid and are accurate 
since e is relatively small. 

5.4.2 Power-law distribution 

For the power-law distribution 



we obtain 



2 a + 2 

X ~ V 2 . 

a + 2 



fo ~ v 



Vf{t) 



-a(a + 2)t 



a + 2 



(103) 



(104) 



and 



9(t) ~ - s , 

implying that the criterion ()50|) is fulfilled. The front profile is given by 



(105) 



u(v,t) ~ $ 



>/2- 



q + 2 
f 2 



J\a{a + 2)t 



(« + 2)tV2 

Concerning the validity of this approach in relation with the criterion ()49j1 . we find that e 



(106) 



J(a + 4) 



Itt 



3a(a+2) 



for t — > +CX3 so that our approximations are marginally valid. 



6 A class of exactly solvable Fokker-Planck equations 

In this section we give a class of Fokker-Planck equations (|18|) for which our approach turns 
out to be exact. More precisely, we derive a relationship between D(v) and U(v) under which 
the corresponding Fokker-Planck equation (JTSj) can be exactly solved. First, we recall that 
equations (fTHjl and ()32j) are equivalent and no approximation has been made in the derivation 
of Eq. (}3*2*|) from Eq. (|18p. However to pass from Eq. ()3*2*|) to Eq. ()3*3*j) for a general velocity 
field V(v), we have made an approximation and have assumed that the concerned velocities are 
close enough to the position of the front Vf(t). This is in fact the only approximation in the 
approach developed above. We now look for a situation where such an approximation is exact. 
This is the case if V(v) = Ax + B where A and B are arbitrary constants. Using 



^(V(v)) = V'(v)^W) = A, 
and replacing V by its expression ()23j) . we equivalently get 



£/<(„)- d —L -l (1„ B)<(„) 
v 2 



A. 



(107) 



(108) 
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Defining R(v) = y/D(v), we obtain 



dv 



U'(v) 



d- 1 



R - R'(v) 



A. 



(109) 



If U(v) and R(v) = y/D(v) satisfy Eq. f)109j) . as discussed in Appendix |D1 we can use the 
results of Appendix [X] with g(t) = A to obtain exact solutions of Eq. (jl8j) for all times. 
Let us make these solutions more explicit. First of all, for h{r) = A, we have 



H(t) = A(t-1), H x {r) = e 



;iio) 



X 2 (r) = \{e A ^ - 1), H 2 (r) = 1(1 - e"^ 1 )). 

We note that the criterion (jSOj) is met only for A < 0. On the other hand, integrating 

dx 



All) 



dt 



V(v) = Ax + B, 



with Xf(l/2) = 0, we get 



x f (t) = |( e 4<*-i> _ 1). 



;ii2) 



(113) 



Then, substituting in Eq. ()152j) . we obtain the general solution 



u(v, t) 



-2At\ 



-00 {e~ At {x- | (e A '-l)]-!/> 2 



u(y,0)dy, 



where we have taken the origin of times at t = 0. If we start from a step function u(x, 0) 
the above expression reduces to 



u(v,t) = $ 



x- f (e A *- 1) 

2 ( e 2At _ ^1/2 



On the other hand, for A = 0, we get 

/i(r)=0, if(r) = 0, #i(r) = l, # 2 (r) = r-l, 



(114) 

T](x), 

(115) 



;ii6) 



The general solution is 



X 2 (r) = 2(r-1), x / (t) = |(2t-l). 



1 r+°° lx - B t-v)* 
V4vrt J-oo 



and if we start from a step function u(x, 0) = r](x), we get 



(117) 



1181 



;ii9) 
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We recall that the above profiles are the exact solutions of the Fokker-Planck equation (|18j) 
when the functions U(v) and D(v) satisfy the differential equation (|109J) . Let us examine some 
particular cases. If A = 0, Eq. ()1(J9J) can be integrated at once leading to the relation 

D{v)=B 2 U w d ~ l e- u(w) dw\ ^zr y (120) 

For the wide class of potentials U (v) satisfying the following condition: 

( d _ i)EM _ u'Hv) < U'{v) 2 (121) 
v 

for large v, the behavior of the diffusion coefficient at infinity is given by 

In particular, when U(v) = v 2 /2 the diffusion coefficient D(v) behaves like C/v 2 for v — > +oo 
and when U(v) = 7t> it tends to a constant. Note that condition (j!21j) is not only satisfied 
by any non constant polynomial, but also by potentials U(v) that dominate ln(t>) for large 
velocities, as for instance U(v) = (\n(v)) a , with a > 1. Note finally that for v — > 0, the 
diffusion coefficient behaves like l/v 2 ^ -1 ' which is divergent except for d = 1. 

In the case of a quadratic potential U(v) = v 2 /2, the integral in Eq. ()120|) can be expressed 
in terms of the error function. There is a simplification in d = 2 leading to D(v) = 1/v 2 . This 
is precisely the case that we have found in Sec. 15.1.11 In that case, we have x = v 2 /2 so that 
Eq. (1118(1 can be explicitly written in terms of v. Inversely, for D(v) = 1/v 2 and g = A, we 
find from Eq. (ITUH1) that U(v) = (A/8> 4 + (B/2)v 2 + (d - 2)lnv so for this type of potential 
the solution of (| 11 4(1 applies with x = v 2 /2. In the case of a linear potential U(v) = ^yv, the 
integral in Eq. (jl20() can be expressed in terms of the incomplete Gamma function. There is 
a simplification in d = 1 leading to D(v) = 1. In that case, we have x = v so that Eq. (|118j) 
can be explicitly written in terms of v. Inversely, for D(v) = 1 and g = A, we find from Eq. 
flEm that U (v) = (A/2)v 2 + Bv + (d-l)\nv so for this type of potential the solution of (ITT^ 
applies with x = v. 

For d = 1, D(v) = 1 and U(v) = (A/2)v 2 + Bv, we have V(v) = Ax + B with x = v 
and we are in the situation mentioned above. If we take f(v,0) = 8{v — vq) and recall that 
u(v,t) = f(v,t)/f e (v), Eq. ()114j) gives after simplification 



1 



f(v,t) = e -w=^«~»*- M +* x -*- M W (123) 

In particular, for 5 = and A = 1, we recover the well-known solution PQ: 

f(v,t) = e n^T. (124) 



e 



— 2*> 



In that case, (v)(t) = v e 1 . Alternatively, for A = and B — 1, we get 

1 (ti-Un+i) 2 

/( V ,t) = _= e (125) 

V47rt 

In that case, (v)(t) = vq — t. Note that there is no normalizable stationary state for the 
Fokker-Planck equation (|TK(l when U(v) = Bv so that f(v,t) tends to zero for large times and 
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spreads so as to conserve the total mass. By contrast, when A ^ 0, the distribution function 
relaxes towards f e (v) for t — > +00. From Eq. ()123|1 . we can obtain the time evolution of the 
distribution function for any initial condition fo(v) = f(v, 0) by multiplying Eq. (|123|) by fo{vo) 
and integrating over v . 

Note finally that for d > 1, the velocity v = |v| is restricted to positive values so that the 
preceding results must be slightly revised. The idea is to extend f(v,t) by parity to negative 
values of v. We replace D(v), U(v) and v d ~ l in Eq. (fTHj) by D(\v\), U(\v\) and \v\ d ~ l . With 
this extension, Eq. (JTSJl is invariant by the transformation v — > —v. Therefore, if f(v,to) 
is even, f(v,t) will remain even for all times and its values for v > correspond to those 
of the distribution function that is solution of the original Eq. (j!8j) . Thus, we can apply 
the preceding results with almost no modification. We note however that since V(v) and 
x(v) = Jq dwj \jD{w) are odd, the case V(v) = Ax + B is only possible if B = 0. By contrast, 
in d = 1, we can consider cases where f(v,t), D(v) and U(v) have no special parity. 

7 Conclusion 

In this paper, we have developed a general formalism to characterize the evolution of the 
distribution function tail for systems described by a Fokker-Planck equation with a diffusion 
coefficient and a friction force depending on the velocity. Our analytical results give good 
agreement with the numerics even in cases where the validity of our approach is marginal. 
When the diffusion coefficient decreases algebraically with the velocity, the progression of the 
front is also algebraic. When the diffusion coefficient decreases like an exponential or like a 
Gaussian, the progression of the front is logarithmic. The high velocity component of the 
distribution function keeps the memory of the initial condition for a long time and is slowly 
depleted. There are several applications of our formalism to, e.g., stellar dynamics, plasma 
physics, vortex dynamics, the HMF model, optical lattices etc. In future works, we shall study 
more specifically the dynamics of point vortices and investigate the evolution of the front profile 
and the time evolution of the correlation functions [S]. 
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A General solution of Eq. (1351) 



In this Appendix, we provide the general solution of the PDE: 

36 d 2 6 , „ . 86 , 

for an arbitrary initial condition (f>i(z) = 6(z, 1). Taking the Fourier transform of Eq. ()126|) 
with the conventions 

6{z) = I HOe^dC, 0(0 = / M*-* 1 ^ (127) 



and using the relation 



we get 



^ = (Hr)-e^ + h(r)^. (129) 

We introduce the change of variables 

f(y,T) = 4>(H 1 (r)y,r), i = H x {r)y (130) 
and choose the function H x {t) such that 

ffi(r) " 2 • [U1) 
Substituting Eq. (fTSTIjl in Eq. (fP23j) . we find that /(j/,t) satisfies 

^ + ^iW-^)]/ = 0. (132) 
Let H(t) be the primitive of h(r) such that 

#(r) = J h(r')dr'. (133) 

Then, we choose H 1 , solution of Eq. (|131|) . such that 

H 1 (r)= e - ^. (134) 
By convention, H(l) = and -Hi(l) = 1. Equation ()132j) can be integrated leading to 

f(y,r) = f(y, \)e^ e^ H ^\ (135) 



where we have defined 



H 2 (t) = J* H x {r') 2 dr' = J* e'^dr'. (136) 
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Returning to original variables, we obtain 



4>{t, r) = 0! [jX^j e^e~^ e . (137) 



We now observe that 

HKr) 



f\ H ^~ H ^dr' = \x 2 (t), (138) 



where x( r ) is t ne function introduced in Eq. (J39|) . Therefore the general solution of Eq. (|126j) 
in Fourier space can be written 



Defining 



where 



q(z)=MH 1 (r)z) <-> = ] (140) 



g(z) = G{z/ X {t)) <-> g(0 = X (t)G(x(t)0 (141) 



G{z) = e-* 2 <-> G (£) = J_ e -? 2 /4 (142) 



we can rewrite Eq. (|139|) in the form 



kM = ^mm- (143) 



Taking the inverse Fourier transform, we can express the solution of Eq. ()126|) as a convolution 

2^ 



or, equivalent ly 

-+oo 



f dz' 
J q(z~z')g(z') — , (144) 



1 P+OO 

z,t) = -= e- x2 MHi(r)(z - X {r)x))dx. (145) 

V 71 " J-oo 



By direct substitution, we can check that Eq. ()145j) is indeed solution of Eq. 

If 4>i(z) = rj(z) is a step function with r)(z) = 1 for z < and 77(2) = for z > 0, we 
immediately find that 

r) = 4= e^'dx = $ (-^) , (146) 

and we recover the result of Sec. El Then, the general solution of Eq. ()126j) can be put in the 
form 



z \ 1 



"+00 

^-•") = (1> ( ^yj y e-^^-^^^^-x^x))^. (147) 
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We introduce the new variable z' = z/x{t) and consider the limit t — ► +00. The integral 
depending on the initial condition is given by 

1 r+oc 

J = -p/ e- I ' 2 (0 1 -r / )(v / 2H2F)(^ -x))dx. (148) 

V 71 " J-oo 



We shall assume that 4>i(z) — > 1 for z — > —00 and <pi(z) — > for 2 — > +00. Then, if #2(7" 
+00 for t — > +00, the function (0i — tj)(^2H2(t)(z' — x)) will be very peaked around 
and we can approximate the integral by 



x = z' 



g-2' 2 /H-'X. 



n 



r])(x)dx-^0. (149) 



Therefore, the condition that the solution of ()126|) tends asymptotically to the function ()146|) 
for r — > +00 is that (0i — rj) (±00) = and 

H 2 (r) = /V^'V -> +00 (150) 



for r — > +00. We can also write the general solution ()145|) in the form 

1 



+00 

z.t) = —= = I e~ ' ' 6Jx)dx. (151) 



Returning to original variables, we get 



1 r+°° (H 1 (2t)( 3: - 3:/ (t))- H ) 2 

u(v,t) = / e 2 «2(2*) u(y,l/2)dy. (152) 



^/27iH 2 {2t) J- 

where a; = x(t>) and we have taken by convention Xf(t = 1/2) = 0. 



B Asymptotic behavior of Eq. (1631) 



Setting z = (7 /2)y 2 , we can rewrite Eq. (|63|) in the equivalent form 



1 7.2 ft* e g -% v * 



X ~ 7y^ e ^' / — ( 153 ) 



With the change of variables y = — z + ?f 2 , we get 



1 7.2 f^ 2 e"" 



x e * v / riy. (154) 



fV 2 



Then taking the limit u — > +00, we find that 



1 /-+oo 1 

J- Z„2 / _„, , 1 7„2 



x e *« / e ~^y ei v . (155) 

7^ Jo 'jv 
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C Criterion for the monotonicity of f(v,t) 



In this Appendix, we establish a criterion which guarantees the monotonicity of f(v,t) for all 
times, if the distribution function is initially monotonic (decreasing). Let us first rewrite the 
Fokker-Planck equation (|18jl in the form 



df ^d 2 f 
— = D — - + 

dt dv 2 



D' + DU' + 



d 



av 



{DU')' + d —^DU' 



(156) 



Because of the positivity of the diffusion coefficient, a classical comparison principle states that 
if f(v,t) is initially positive, it will remain positive for all times [23J. The idea is now to apply 
the same argument to g = df /dv. Taking the derivative of Eq. (j!56|) with respect to v, we get 



dg_ 
dt 



If 



then 
dg 



D 



8Fg_ 
dv 2 



+ \2D' + DU' 



d- 1 



D 



dg_ 

dv 



+ 



D" + 2(DUJ + - — -DU' + (d-l) 



+ 



(DU')" + (d-l) 



DU' 



(DU')" + (d-l) 



DU' 



<0, 



/■ (157) 



(158) 



<D*l + 
dt dv 2 



2D' + DU' 



d-l 



D 



dg_ 

dv 



D" + 2(DU')' + 



d 



—DU' + {d-l) (- 
v V v 



9- 

(159) 



Again we use a comparison principle and deduce from this inequality that if g(v,t) is initially 
negative, it will remain negative for all times. Therefore, if f(v,t) is initially decreasing, it 
will remain decreasing for all times if the criterion ()158j) is satisfied. This criterion is just a 
sufficient condition of monotonicity. 

Let us give some particular examples of application. In the case of a quadratic potential 



U(v) 



2 /2, the criterion (|158J) becomes 



(Dv)" + (d - 1)D' < 0. 



(160) 



For a diffusion coefficient of the form D = v~ a , we find the condition a < d. In particular, 
for the case a = d = 3, the criterion is satisfied. We can check that the criterion ()16U|) 
is also satisfied for all v for the diffusion coefficient (JJJ-a) corresponding to the Coulombian 
(or Newtonian) potential of interaction. Alternatively, for a diffusion coefficient of the form 
D = e~ lv , we find that the criterion (J16U)) is not satisfied for large values of v (more precisely 
v > (d+ a/g? 2 + 87)747). Therefore, the profile f(v,t) can be non-monotonic even if the initial 
condition is monotonic, as in the case of Fig. El 



D General solution of Eq. (11091) 



In this Appendix, we show that Eq. (|109J) can be solved explicitly. First of all, we note that it 
can be written 

U'(v)-^-—]r-R'(v)=Ax + B. (161) 
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Since R— dv/dx, we obtain 



U'(v) 



d-1 R'(v, 



v R(v) 



dv = (Ax + B)dx 



which leads to 



U{v)d-l 



e v 'v 



R{v) = Ke~^ x+Bx) 



where if is a constant. Equation (j!63|) can again be integrated into 



+00 



e -u {w ) w d-i dw = K e-^ y2+B ^dy. 

J x(v) 



For A = 0, we obtain 



(162) 



(163) 



(164) 



+00 



- u ^w d ~ l dw 



K 



-Bx 



(165) 



Substituting this relation in Eq. ()163|) . we recover the result of Eq. (jl20j) . On the other hand, 
for A ^ 0, we get 



V A 



<A( B 
2 I + I 



(166) 



where $(x) is defined in Eq. (jHEJ)- Substituting the foregoing relation in Eq. ()163|) . we find 
that 



D(v) 



e 2U(v) 
v ^d-l) ■ 



Ce~ u{w) w d - l dw 



where 



F(x) = exp{-2($- 1 (a;)) 2 } 



(167) 



168) 



and C is a constant. 
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